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We studied the (00|) diffraction peak in the low-temperature phase of magnetite (FeaCU) us- 
ing resonant soft x-ray diffraction (RSXD) at the Fe-Z/2,3 and O-K resonance. We studied both 
molecular-beam-epitaxy (MBE) grown thin films and in-situ cleaved single crystals. From the com- 
parison we have been able to determine quantitatively the contribution of intrinsic absorption effects, 
thereby arriving at a consistent result for the (00 1) diffraction peak spectrum. Our data also allow 
for the identification of extrinsic effects, e.g. for a detailed modeling of the spectra in case a "dead" 
surface layer is present that is only absorbing photons but does not contribute to the scattering 
signal. 

PACS numbers: 71.30.+h, 61.05.cp 



Magnetite, Fe3C>4, shows a first order anomaly in the 
temperature dependence of the electrical conductivity at 
TV ~ 120 K, the famous Verwey transition!^ It is accom- 
panied by a structural phase transition from the cubic 
inverse spinel to a distorted structure. One usually con- 
nects this transition with charge ordering of the Fe 2+ and 
Fe 3+ ions on the octahedrally coordinated, so-called B- 
sites. In a recent diffraction study, Wright, Attfield and 
Radaelli found long-range charge order from the pattern 
of shorter and longer bond lengths between B-site iron 
and oxygen ions below TypEl Subsequently, using this 
structure as input, local density approximation + Hub- 
bard U (LDA+U) band structure studies have calculated 
the corresponding orbital ordering, which involves mainly 
the minority spin electron in the t^g state of the 2+ B- 
sitesPH 

Recently this ti g orbital order has been studied using 
RSXD at the O Is -> 2p (K) 6 and Fe 2p 3d (L 2 , 3 fi 
resonance. The O-K edge resonance enhancement of the 
(00 1) diffraction peak (notation refers to the cubic room- 
temperature unit cell) on an ex-situ polished bulk single 
crystal was interpreted by Huang et al. as a signature 
of a particular charge/orbital order at the oxygen-sites P 
At the Fe L2,3 _ edges, using MBE-grown thin films, the 
(00 1 ) maximum energy coincides with the resonance of 
the 2+ B-site ions and was assigned by Schlappa et al. 
to Fe ^g-orbital order on these sitesP 

Very recently these results have been challenged. 
Wilkins and coauthors have presented RSXD data from 
the (00|) peak [(001) in the orthorhombic notation used 
in Ref. [8] of bulk magnetite at the O-K and Fe-L2,3 
resonanceP. Their O-K spectrum is similar to that of 
Huang et alW but their Fe-L2,3 data are very differ- 
ent from those of Schlappa et aZP Almost at the same 
time Garcia et a^P reported a Fe-L 2 ,3 RSXD spectrum 
from a bulk sample, which is very different from those of 



both Wilki ns e t al. and Schlappa et al.. In both recent 
publications^! the shape of the respective Fe-L2,3 spectra 
are explained by absorption effects. Presently, hence, at 
least three different Fe-L 2j 3 RSXD spectra for the (00^) 
peak of magnetite have been published and very different 
conclusions have been drawn about the electronic origin 
of the peak. 

Here we report our efforts to resolve the confusion 
about the RXSD data on magnetite. Our strategy is 
to carry out the experiments using a magnetite bulk sin- 
gle crystal which was carefully cleaved in-situ in order to 
obtain the best possible surface quality, and to compare 
the results with the measurements on high quality MBE- 
grown magnetite thin films. In this manner we are able 
to distinguish intrinsic absorption effects from extrinsic 
effects. Taking intrinsic effects properly into account, we 
find consistent results for in-situ cleaved bulk samples 
and thin films. Furthermore, taking a so-called "dead" 
surface into account, which leads to extrinsic absorption, 
we are able to provide a quantitative explanation for the 
varying RSXD spectra reported in the literature so far. 

The soft x-ray scattering experiments were performed 
at the HZB beam line UE46-1 at BESSY using the two- 
circle UHV diffractometer designed at the Freie Univer- 
sitat Berlin. We used 40 nm thick magnetite films epitax- 
ially grown on an epi-polished MgO(OOl) substrates and 
several magnetite bulk single crystals. The samples were 
oriented such that two base vectors of the cubic room- 
temperature unit cell (a = 8.396 A) were parallel to the 
scattering plane. The polarization of incoming light was 
either parallel or perpendicular to the scattering plane 
(tt- or cr-polarization, respectively). A silicon-diode pho- 
ton detector without polarization analysis was used with 
the angular acceptance set to 1° in the scattering plane 
and 5° perpendicular to it. The energy dependence of 
the diffraction-peak intensities (diffraction spectra) from 
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FIG. 1: (color online) RSXD spectra of the (OOf) peak of 
bulk magnetite (open symbols) and magnetite thin film (filled 
symbols) at (a) Fe Z/2,3-edges and at (b) O if -edge. The 
spectra are scaled to equal height. 



FIG. 2: (color online) RSXD spectra of the (00±) peak of 
bulk magnetite without (open symbols) and with absorption 
correction (filled symbols) at (a) Fe 1,2,3-edges and at (b) O 
K-edge. The spectra are scaled to equal height. XAS spectra 
are shown as solid lines. 



the film sample were recorded by varying the photon en- 
ergy and keeping the momentum transfer constant. For 
bulk samples this method leads to distorted spectra be- 
cause of the much narrower peaks and refraction effects. 
We therefore performed longitudinal (L) scans at every 
photon energy and integrated the (OO^)-peak intensity. 

Fig. Ill displays the (00 \) resonance spectra at the Fe 
£2,3 and at O if-edges for bulk magnetite and magnetite 
thin film. One can observe very distinct differences be- 
tween the bulk and the thin film. At the Fe edge, for 
example, the first peak at 707 eV is much more pro- 
nounced in the bulk than in the thin film. At the O 
edge, the energy position of the peak of the bulk is clearly 
shifted with respect to that of the thin film. We ascribe 
these differences to a change of the scattering volume due 
to variations of the penetration depth of the light as a 
function of its energy. This effect is quite severe for en- 
ergy variations across resonances in the soft x-ray regime 
due to the very high photo-absorption coefficients. For 
a bulk sample, the scattered intensity in specular geom- 
etry is proportional to with fi being the absorption 
coefficient.^ This absorption effect can be accounted for 
by multiplying the resonance spectra by fi; fi can be re- 
liably obtained from the x-ray absorption spectroscopy 
(XAS) signal recorded in the total electron yield mode. 

The effect of this correction is shown for bulk mag- 
netite in Fig.[2]for both the Fe-L2,3 and O-K resonances. 
The open symbols are the raw data and the filled sym- 
bols are those after correction. The absorption signals are 
shown as solid lines. The correction for the Fe-spectrum 
affects mainly the relative intensities of the spectral fea- 



tures, while in the O-spectrum, the main peak, which is 
just on the slope of the absorption signal, becomes shifted 
in energy. 

For a film of thickness D the factor to correct for ab- 
sorption effects is fi/(l — exp(— 2 fid)), d — D/ sin9 de- 
scribes the effective photon path length at incident angle 
0. To perform this correction requires to know fi in ab- 
solute units. Because of the thinness of the films, the 
correction should be smaller than for the bulk. 

In Fig. [3] we now compare again the corrected spec- 
tra from the bulk (filled symbols) with the uncorrected 
film spectra as presented in Ref. [7] (open symbols) . For 
the O-K spectra the agreement between bulk and film 
data is now close to perfect. Also for the Fe resonance 
the overall agreement between the spectra is very good, 
all spectral features agree. The main discrepancy here is 
in the relative intensities of the L3 and L2 part of the 
Fe-resonance spectra. The application of an absorption 
correction on also the thin film data would probably re- 
solve much of this discrepancy. Thus, apart for this, we 
find very similar spectra from the film and the bulk at 
both resonances. We can, hence, safely conclude that 
the resonant diffraction data recorded from the thin film 
sample are representative for the material magnetite as 
such and that the conclusions derived from the data of 
magnetite thin films in our former Letter^ stand for the 
intrinsic properties of magnetite. 

We now compare our experimental data with those 
from Huang et alP, Wilkins et alW, and Garcia et alW 
from bulk magnetite. We first focus on the Fe i^^-edges 
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FIG. 3: (color online) RSXD spectra of the (00 1) peak from 
magnetite thin film and bulk crystal at (a) Fe Z/2,3-edges and 
(b) O K-edge using n polarized light. 
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FIG. 4: (color online) Fe Z/2,3-edges RSXD spectra of the 
(00 1 ) peak and XAS spectra taken from our magnetite bulk 
crystal (red) and from Wilkins et aZ.P (black). The green 
thick line is the result of a simulation multiplied by a fac- 
tor 500 (see text). The XAS spectra are included as energy 
reference. 



data. In Fig.[4]we present our data and those reproduced 
from the paper of Wilkins et aZ.P The energy scales were 
aligned using the XAS data to account for small energy- 
scale offsets between different soft x-ray monochroma- 
tors. The energy scale given at the top of Fig. prefers to 
the one used by Wilkins et aiP, and the bottom scale is 
the one of our experiment. It was claimed by Wilkins et 
al. that their data were strongly distorted by absorption 
effects, which suppressed all scattered intensity when the 
absorption increases. However, this cannot be the case 



because our raw data are already very different. In fact, 
it is quite unlikely that absorption in an otherwise intact 
sample can suppress the scattered intensity completely.^ 
What could cause such strong distortions, however, is 
a layer on top of the sample, which does absorb photons 
but does not contribute to the scattered signal. This 
could be either a polycrystalline layer or a crystalline 
layer with different surface orientation as the bulk or a 
layer of different chemical composition. That this " dead" 
layer scenario is likely is demonstrated in Fig. [4j The dia- 
mond symbols are the experimental data from Wilkins et 
alPand the green thick solid line is a simulation in which 
we describe the presence of an absorbing layer of thick- 
ness D a by multiplying our experimental RSXD spec- 
trum by e -2d al j.(hv) with da = £> o / s i n 0. 

We find the best agreement between data and simula- 
tion for d a X ji m 18 at the absorption maximum. The 
resulting spectrum reproduces not only the peak posi- 
tion but also the shape of the data from Wilkins et aZP 
very well. We note that as compared to our experimental 
spectrum the green line has been multiplied by a factor 
500; more than 99 percent of the intensity of the diffrac- 
tion spectrum is taken away by the absorbing dead layer. 
Concerning the Fe L 2j 3 data from Garcia et oPwe note 
that their (raw) data are closer to ours, although there 
are appreciable quantitative differences. It therefore ap- 
pears likely that the effect of an absorbing dead layer is 
present also in their data, although of a much thinner 
one than in the work of Wilkins et alW 

We now address the oxygen resonance spectra. Experi- 
mentally both in the papers by Wilkins et aZP and Huang 
et alPa single resonance peak is observed, while our data 
show clearly a satellite 1.3 eV higher in energy. Huang 
et al. observe additionally a broad hump at about 4 eV 
higher photon energy, which is missing in our data. The 
1.3-eV satellite occurs at energies where the absorption 
cross section is high and the filter effect discussed above 
for the Fe resonance may suppress this feature. Since no 
absorption data for energy reference are given by Wilkins 
et a/P we focus on those from Huang et aZ.P. In Fig. [5] 
(a) we reproduce the O-K RSXD data from Huang et 
alW. Again, we aligned the energy scales using the O-K 
XAS spectra (thin/dotted lines). Since our experiment 
was performed with higher energy resolution, we broad- 
ened our experimental XAS data slightly to obtain the 
same energy resolution as in the experiment of Huang et 
alP. On such a calibrated energy scale we find the res- 
onance peak at a by 230 meV higher photon energy, i.e. 
much closer to the absorption maximum than Huang et 
al. did. 

We have carried out the same simulation as in the Fe 
^2.3-edges analysis above for O iC-resonance as shown in 
Fig. [5] (a). The triangle symbols are the experimental 
data from Huang et a/P and the green thick solid line is 
the simulation result. The resulting spectrum reproduces 
the main peak position and shape, and also contains the 
weak and broad hump at 4 eV above the main peak, i.e. 
a structure which is not present in the real spectrum and 
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FIG. 5: (color online) (a) O K-edge RSXD spectra of the 
(00 1) peak and XAS spectra taken from our magnetite bulk 
crystal (red) and from Ref. [B] (black). The green thick line is 
a simulated result (see text). The XAS spectra are included 
as energy reference, (b) Red line and symbols is our RSXD 
spectrum, and green thick line is the same simulated result as 
in (a) but now both on the same vertical scale, showing that 
more than 95 percent of signal is taken away by the absorbing 
layer. 

which is artificially generated by the local minimum in 
the absorption at that photon energy. In panel (b) we 
plot our original spectrum (red line and symbols) and 
the one generated by assuming an absorbing layer (green 
solid line) on the same vertical scale. One can also see 



how much intensity of the diffraction spectrum is taken 
away by the absorbing dead layer: more than 95 percent 
of the signal is gone. The shape of the spectrum of Huang 
et alPis completely destroyed, the 1.3-eV satellite is sup- 
pressed, and the 4-eV feature is artificially generated by 
the dead surface layer. 

We note that assuming a photon mean free path (1/m) 
of about 1300 A (a value proposed by Huang et aZP) 
at the energy between 528 and 531 eV, and considering 
the detection geometry with photon incidence and detec- 
tion angle around 45 degrees, we find a thickness of the 
"dead" surface layer of about a quarter of a micrometer. 
A similar thickness comes out for the data of Wilkins et 
aZP at the Fe-L3 when we estimate the photon mean free 
path at the Fe-L3 maximum based on our XAS data. 

In summary, we have recorded RSXD spectra of the 
(00 \ ) peak at both Fe L 2 ,3-edges and O K-edge from in- 
vacuo cleaved magnetite single crystals, which essentially 
agree with those from magnetite thin films published ear- 
lier. We found that our spectra differ strongly from those 
published by Huang et alP and Wilkins et aZP, and to 
lesser extent, from Garcia et aZP We were able to as- 
cribe the distortions in their spectra to extrinsic effects, 
namely the presence of a dead surface layer. Both our 
bulk and thin film data can hence be used for a reliable 
modeling of the charge and orbital ordering phenomena 
in magnetite. 
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P. Abbamonte, private communication: For a single res- 
onance the scattered intensity I = [/' + if"\ 2 at the res- 
onance maximum is proportional to fi 2 , because the real 
part of the scattering amplitude, /' is zero at the resonance 
maximum and /" oc /i. The absorption effect, however, is 
linear in fi, which means that the observed intensity will 
still be finite and proportional to /i. While the actually 
probed linear combination of scattering amplitudes may 
weaken this simple argument, it is still very unlikely that 
absorption can reduce the scattered intensity to zero for a 
wide range of photon energies. 



